Vacuum Rabi oscillations in a macroscopic superconducting qubit LC oscillator system 
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We have observed the coherent exchange of a single energy quantum between a flux qubit and a 
superconducting LC circuit acting as a quantum harmonic oscillator. The exchange of an energy 
quantum is known as the vacuum Rabi oscillations: the qubit is oscillating between the excited state 
and the ground state and the oscillator between the vacuum state and the flrst excited state. We 
have also obtained evidence of level quantization of the LC circuit by observing the change in the 
oscillation frequency when the LC circuit was not initially in the vacuum state. 

PACS numbers: 



The coupling of quantum two-state systems to har- 
monic oscillators constitutes a broad and active research 
field, most notably in atomic cavity quantum electro- 
dynamics (QED) where atoms having a large electrical 
dipole moment interacts with the quantized electromag- 
netic field in high-Q cavitiesi, but also in other systems, 
for example, ions held in linear traps where the ions cou- 
ple to a vibrational mode acting as a harmonic oscillatou^, 
and in solid state, where coupling of atom-Hke objects to 
harmonic modes has been pursued for some time, start- 
ing from quantum wells in micro-cavities^ to recent ex- 
periments on quantum dots in photonic crystals^, and 
two experiments with superconducting devices: a system 
with a fiux qubit coupled to a superconducting quan- 
tum interference device (SQUID) showed Rabi oscilla- 
tions driven by an external fieldS^ and the vacuum Rabi 
splitting due to strong coupling of a Cooper pair box 
to a microwave photon in a superconducting waveguide 
was demonstrated^. These systems attract increasing at- 
tention because of their potential use in future quantum 
information processing. In addition, they offer an ideal 
testing ground for studying the fundamental interactions 
between matter (atoms, qubits) and light (harmonic os- 
cillators) . 

In this letter we report the observation of vacuum Rabi 
oscillations in a macroscopic superconducting solid state 
system. In our strongly coupled system, the qubit state 
and the qubit interaction time with the oscillator are 
controlled by a combination of microwave and DC-shift 
pulses, resulting in a measuring sequence analogous to 
atomic cavity QED. We demonstrate vacuum Rabi oscil- 
lations between two quantum objects with macroscopi- 
cally distinct statesS.: the states are the coherent motion 
of the superconducting condensate involving millions of 
Cooper pairs, or the order of a few 100 nA in terms of 
current associated with the states. 

Superconducting Josephson junctions based qubits are 
ideal artificial atoms where quantum coherence has been 
demonstrated for many different typesS*ifliii*iSii2iiiii&. 
Recently, two qubit operations in the time domain have 
also been demonstratediS*ii»i^. In our experiment a su- 
perconducting persistent current qubit plays the role of 




Figure 1: (a) SEM micrograph of the sample. The qubit and 
the detector SQUID enclosing it are the small square loops in 
the lower center picture. The square plates at the top of the 
picture are the top plates of the on-chip capacitors separated 
by an insulator from the large bottom plate, (b) A close up 
of the qubit and the SQUID. The qubit dimension is 10. 2x 
10.4 iJLva? . (c) Equivalent circuit of the sample. The .losephson 
junctions are indicated by crosses: three in the inner qubit 
loop and two in the SQUID. The LC mode is indicated by the 
dashed line. The inductance and capacitances are calculated 
from the geometry to be L = 140 pH and C = 10 pF, and 
the qubit LC oscillator mutual inductance to be M = 5.7 pH. 
The current and voltage lines are flltered through a series 
combination of copper powder fllters and lossy coaxial cables 
at mixing chamber temperature and on-chip resistors [Rv = 
3 kfi and Ri = 1 m). 



the atom. The persistent current qubit is a superconduct- 
ing loop interrupted by three Josephson junctions (see 
Fig. 1) with one junction smaller by a factor of a = 0.8. 
In the vicinity of half a fiux quantum in the loop the de- 
vice is an effective two-level systemic described by the 
Hamiltonian i?qb = h/2 {tOz + Actj:), given in the basis 
of clockwise and counterclockwise currents. Here <Jz/x 
are the Pauli spin matrices, he = 2Jp ($ex — $o/2) is the 
energy bias (/p is the persistent current in the qubit, <i>cx 
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Figure 2: (a) Spectroscopic characterization of the qubit- 
oscillator system showing the LC oscillator at Vv = 4.35 GHz 
and the qubit dispersion around the gap of A = 2.1 GHz. (b) 
A close up of the region around 4.35 GHz (indicated by the 
red square in (a) showing an avoided crossing. The lines are 
guides for the eye. (c) Schematic of the pulse sequence used 
to obtain the spectroscopy in (b): the qubit operating point is 
fixed at 10.5 GHz and via the MW line a shift pulse of variable 
height moves the operating point to the vicinity of 4.35 GHz. 
Two 50 ns MW pulses separated by 2 ns are added to the shift 
pulse. Here the second MW pulse is phase shifted by 180°. 
The phase shifted second pulse damps the oscillations in the 
LC circuit-fi, and is crucial in terms of resolving the relatively 
weak qubit signal in this region. After the MW pulses the 
qubit state is measured by applying a measurement pulse to 
the SQUID (green curve). The spectroscopy in Fig. 2(a) was 
obtained with the same scheme, but without the second phase 
shifted pulse. 



is the external flux threading the qubit loop, and $o is 
the flux quantum), and A is the tunnel splitting. The 
energy gap of the q ubit, controlled by the external flux, 
is hE = We2 + A2. 

In our sample geometry, the qubit is spatially sepa- 
rated from the rest of the circuitry. The qubit is en- 
closed by a superconducting quantum interference device 
(SQUID) that is inductively coupled to the qubit. The 
SQUID functions as a detector for the qubit state: the 
switching current of the SQUID is sensitive to the flux 
produced by the current in the qubit. The qubit is also 
enclosed by a larger loop containing on-chip capacitors 
that provide a well-deflned electromagnetic environment 
for the SQUID and flltering of the measurement leads. 
The lead inductance L and capacitance C in the outer 
loop constitute an LC oscillator [see Fig.^b)] with reso- 
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Figure 3: (a) Vacuum Rabi oscillations (symbols) and a nu- 
merical fit (solid line), (b) The few lowest unperturbed and 
dressed energy levels when the system is in resonance, (c) 
The qubit energy level diagram and pulse sequence for the 
vacuum Rabi measurements. The tt pulse (4.6 ns long at -16 
dBm) on the qubit brings the system from state 1 to 2 and 
the shift pulse changes the flux in the qubit by $ahift, which, 
in turn, changes the operating point from 2 to 3 where the 
system undergoes free evolution between |e, 0) and |g, 1) at 
the vacuum Rabi frequency SIr until the shift pulse ends and 
the system returns to the initial operating point where the 
state is measured to be either in 2 or 4. 



nance frequency cUr = 27rj^i. = 1/y/LC. The LC oscillator 
is described by a simple harmonic oscillator Hamiltonian: 
-f^osc = (a^a + 1/2), where (a) is the plasmon cre- 
ation (annihilation) operator. The qubit is coupled to 
the LC oscillator via the mutual inductance M, giving 
an interaction Hamiltonian Hi = hXaz (a^ + a) , where 
the coupling constant is hX — MIp^/hLu^j2L. The total 
system is thus described by a Jaynes-Cummingg- type of 
Hamiltonian H — h/2 {ecjz + ^(Jx) + hLOr (a^a -t- 1/2) -|- 
hXa^ {a) -f a) . We denote the state of the system by 
\Q,i), with the qubit either in the ground {Q = g) or 
excited {Q = e) state, and the oscillator in the Fock 
state {i = 0, 1, 2, . . .). The parameters of the system can 
readily be engineered during fabrication; the qubit gap 
is determined by a and the junction resistance, the os- 
cillator plasma frequency is fully determined by L and 
C, and the coupHng between the qubit and the oscillator 
can be tuned by M and L. 

The measurements were performed in a dilution re- 
frigerator (T « 20 mK) using strongly flltered mea- 
surement lines and low noise electronics. The system 
was initially allowed to relax to its ground state, where 
the thermal occupation of the excited states is negligi- 
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ble {kBT <C hE,huj^). The readout of the qubit state is 
performed by applying a pulse sequence to the SQUID 
[see Fig. 2(c)], and recording whether the SQUID had 
switched to a finite voltage state or had remained in the 
zero voltage state. The height and length of the pulse 
were adjusted to give the best discrimination between the 
ground and excited states. In our configuration the mea- 
sured switching probability is directly related to the oc- 
cupation of the excited state. Coherent control of the sys- 
tem is achieved by applying microwave (MW) to the sys- 
tem via the MW line (see Fig. 1). The MW line is induc- 
tively coupled to the system, and coherent MW excita- 
tions introduce an oscillating driving term in the Hamil- 
tonian: Hd — Bqh<Jz coswexi + Bosc (a^ + a) coswex^, 
where Bqb (-Bosc) is the amplitude of the qubit (LC os- 
cillator) drive and v^^ = Wcx/Stt is the drive frequency. 
The amplitudes of the drive are proportional to the am- 
plitude of the MW signal, and because of the geometry 
the LC oscillator experiences a four times larger ampli- 
tude, Bqb = 1/4 Bosc. 

We start with a spectroscopic characterization of the 
system to map out the qubit and LC oscillator parame- 
ters. The large scale sweep in Fig. 2(a) shows the qubit 
dispersion curve and LC oscillator at Vr = 4.35 GHz. 
The qubit parameters A — 2.1 GHz and /p — 350 nA are 
extracted from the dispersion curve. Figure 2(b) shows 
spectroscopic measurements obtained in the vicinity of 
the resonant frequency of the LC oscillator (area enclosed 
by the red square in Fig. 2(a) where an avoided crossing 
can be seen. Here, the LC oscillator is strongly excited 
to a classically oscillating state and we cannot draw any 
quantitative conclusion about the coupHng from the spec- 
troscopy, i.e., we are not able to observe the vacuum Rabi 
splitting which usually refers to the splitting of the states 
involving the two lowest oscillator Fock states^. 

Next we investigate the dynamics of the coupled sys- 
tem in the time domain. We performed a measurement 
cycle where we first excited the qubit and then brought 
the qubit and the oscillator into resonance where the ex- 
change of a single energy quantum between the qubit and 
oscillator manifests itself as the vacuum Rabi oscillation 
|e,0) ^ |gi 1) (see Fig. 3). Figure. 3(c) is a schematic 
of the pulse sequence: We started by fixing the qubit 
operating point far from the resonance point [point 3 in 
Fig. IHIc)] and prepared the qubit in the excited state 
by employing a tt pulse. The tt pulse was followed by a 
shift pulse, which brought the qubit into resonance with 
the oscillator for the duration of the shift pulse. After 
the shift pulse the qubit and the oscillator were brought 
back into off-resonance and the measurement pulse was 
applied to detect the state of the qubit. It is important 
to note that the rise time of the shift pulse, Trisc = 0.8 ns, 
is adiabatic with respect to both the qubit and the os- 
cillator, Trisc > [^1 E,2t: / ujj\, but non-adiabatic with re- 
spect to the coupling of the two systems, Tj-iso < 1/A. 
Hence, when the system reaches the resonant point, it 
is in the state |e,0), which is not an eigenstate of the 
total Hamiltonian and therefore free evolution between 



the states |e,0) and |g, 1) begins. The physics behind 
the vacuum Rabi oscillations is thus different from that 
of normal Rabi oscillations where the system is driven 
by an external classical field and oscillates between two 
energy eigenstates. Also, in the normal Rabi oscillations 
the Rabi frequency is determined by the drive ampli- 
tude whereas the vacuum Rabi oscillation frequency is 
determined only by the system's intrinsic parameters. 
The observed Rabi oscillations are in excellent agree- 
ment with those calculated numerically [soHd fine in Fig. 
3(a)]. The numerical calculation uses the total Hamilto- 
nian and incorporates the effects of the decoherence of 
the qubit and the oscillator. The calculation was per- 
formed with the known qubit and LC oscillator parame- 
ters (obtained from spectroscopy and qubit experiments: 
qubit dephasing rate = 0.1 GHz, qubit relaxation rate 
Fe = 0.2 MHz, A = 2.1 GHz, w^/27r = 4.35 GHz) and by 
treating the coupHng constant and oscillator dephasing 
and relaxation rates as fitting parameters. From the fit 
we extracted the coupling constant A = 0.2 GHz, oscil- 
lator dephasing rate F^ — 0.3 GHz, and relaxation rate 
To — 0.02 GHz. The coupHng constant extracted from 
the fit agrees well with that calculated with the mutual 
inductance A = 0.216 GHz. 
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Figure 4: (a) Rabi oscillations when a 2 ns long pulse with 
frequency I'cx = 4.35 GHz and amplitude Amw oc 10^°'^=''^'' is 
inserted between the tt pulse and the shift pulse, (b) Measured 
Rabi oscillations at different drive powers (symbols), and a 
fit (solid curve) to [X]n=o cos{^n + 1 fiat)] exp(— Ft) with 
ao, . . . ,03 and JIr as the only fitting parameters (F is fixed 
from a fit to the -100 dBm curve), (c) The weights of the 
different frequency components ao, • . • ,03 obtained from the 
fit as a function of the drive amplitude (The red arrows show 
the position of the curves in (b). The error-bars indicate the 
errors obtained from the fitting procedure. 
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To confirm that the oscillations we observed originated 
from the coupling of the qubit to a quantized harmonic 
oscillator we performed an experiment where we actively 
excited the oscillator to probe its quantized nature, along 
the lines of Ref. 0. We inserted a short (2 ns) pulse that 
was resonant with the oscillator (4.35 GHz and variable 
power Pose) between the tt pulse and the shift pulse. Fig- 
ure 4(a) shows the Rabi oscillations for increasing power 
to the oscillator. The coherent drive of the oscillator ex- 
cites higher states according to a Poissonian distribution. 
The Rabi frequency involving Fock states n and n + 1 is 
given by \/n + IJIr, hence the Rabi oscillations become 
faster when the oscillator is in higher excited states, as 
seen at the top of Fig. 4(a), and there are several fre- 
quency components in the oscillations, which can be seen 
in Fig. 4(a): a crossover of the vacuum Rabi frequency 
8 : 5.5 w 1 : l/\/2 as the amplitude is increased. In 
Fig. 4(b) the Rabi oscillations are plotted together with 
a fit to an exponentially decaying sum of cosines with fre- 
quency, y^n + lflR and n = 0, . . . , 3. The only fitting pa- 
rameters are the amplitudes of the cosines and the Rabi 
frequency Hr. The decay rate being determined from 
the curve with the lowest power (-100 dBm). The Rabi 
frequency obtained from the fits varied between 137 and 
147 MHz, which is within the accuracy of the fitting. In 
Fig. 4(c) the amplitudes of the cosines are plotted as a 
function of the drive amplitude ^mw oc 
the drive amplitude increases the population of the oscil- 
lator ground state decreases (the populations of the states 
are proportional to the amplitudes of the frequency com- 
ponents) and the excited states become increasingly pop- 
ulated. We observe significant contributions to the Rabi 
frequency from the three lowest oscillator states. With a 
stronger drive to the oscillator the Rabi signal disappears, 
probably due to faster decoherence from the higher ex- 
cited states. In the case of a qubit coupled to a classical 
oscillator there would be one single frequency component 
of the Rabi oscillations, which is not the case in our mea- 
surements where we can map out up to three frequency 
components, showing that our system is in the quantum 
regime in which the discreteness of quantized energy lev- 



els is manifest. When the drive to the oscillator is very 
weak or non-existent the first excited state still has a 
finite population. We attribute this population to the 
strong coupling between the MW-line and the oscillator; 
the oscillator is excited even by off- resonant microwaves. 
This is confirmed by our numerical calculations, in the 
fit of the vacuum Rabi oscillation [Fig. 3(a)] we found 
(0) ~ 0.8 and (1) ~ 0.2 after the tt pulse. Therefore we 
can conclude that the Rabi oscillations observed without 
excitation of the oscillator are predominantly related to 
the ground and first excited state of the oscillator, that 
is, the vacuum Rabi oscillation. With improved design, 
the coupling of the oscillator to the MW-line can easily 
be reduced by more than one order of magnitude. We 
also note that the relaxation and dephasing times were 
shorter for the oscillator than for the qubit. This can 
partly be attributed to the fact that the oscillator was 
galvanically connected to the measurement lines and we 
beheve that in a free standing LC circuit the relaxation 
and dephasing times will be much improved. 

In conclusion ,we have demonstrated vacuum Rabi os- 
cillations in a superconducting soHd-state system, show- 
ing the exchange of a single energy quantum between a 
qubit and a LC circuit. We have also used the qubit 
to probe the state of the LC oscillator and showed that 
the LC circuit displayed characteristics of level quantiza- 
tion. With improved sample design it should be possible 
to generate and measure nonclassical states of the LC 
oscillator, as have been done in ion trap experiments^^. 
Our results supports the suggestions that future quan- 
tum information processing in solid-state devices can be 
achieved with this kind of system, for example coupHng 
several qubits via an oscillato r^^i^Ti^§ . 
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